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We present for the first time numerical results for the particle (photon) creation rate of Dynam- 
ical Casimir effect (DCE) radiation in a resonant cylindrical microwave cavity. Based on recent 
experimental proposals, we model an irradiated semiconducting diaphragm (SCD) using a time de- 
pendent 'plasma sheet' where we show that the number of photons created for the TMon mode is 
considerably enhanced even for low laser powers (of fij order). Conversely to the moving mirror 
case, we also show that the fundamental TM mode (TMoio) is not excited for an irradiated plasma 
sheet. We show that polarization (arising due to the back reaction of pair created photons with the 
plasma SCD) implies losses for TM, but not TE modes. However, we argue that these losses can be 
reduced by lowering the laser power and shortening the relaxation time. The results presented here 
lead support to the idea that TE and, in particular, TM modes are well suited to the detection of 
DCE radiation in a cylindrical cavity. 
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I. INTRODUCTION 

The dynamical Casimir effect (DCE) was first discussed by G. T. Moore [1] more than 40 years ago in 1970, who 
showed that pairs of photons would be created in a Fabry-Perot cavity if one of the ends of the cavity wall moved with 
periodic motion. The number of photons produced during a given number of parametric oscillations is proportional to 
sinh 2 (2wi v/c), e.g., see [2]. However, in most cases, the mechanical properties of the material imply v/c <C 1, where 
v is the wall velocity and c is the speed of light. To overcome this problem there have been various other proposals 
besides the mechanical oscillations of a boundary, such as using a dielectric medium [3-9], also see [10]. This leads 
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FIG. 1. General idea of plasma window experiments such as that in Italy [15-17]: A laser periodically irradiates a semiconductor 
diaphragm (SCD) inside a superconducting cavity. In our simulations we assume a cylindrical cavity with dimensions L z = 100 
mm, R = 25 mm, and we also assume that the SCD can be placed at various locations d within the cavity: rj — d/L z . Inset: 
Typical pulsed laser profile with each pulse train of order 10 ns (10 ~ 100 pulses) repeated every ~ 10 ms. Each pulse is typically 
of order T = t e +to+T~ 100 ps and depends on the excitation time, t e , and relaxation time, t. 

to an effective wall motion by varying the optical path length of the cavity [5, 7, 10]. There are also other methods 
such as illuminated superconducting boundaries [11] and time varied inductance effects in quantum circuit devices 
[12], where possible indirect evidence for photon creation has been reported [13] (also see the very nice review in 
Reference [14]). An experiment in progress [15-17] (see Figure 1 for the general idea) uses a plasma mirror obtained 
by irradiating a semiconductor sheet with a pulsed laser. This leads to an effective wall motion by varying the surface 
conductivity, which generalizes early proposals [18-20] that suggested using a single laser pulse. (For more details on 
all these ideas see the excellent review by Dodonov [21].) 

The goal of this paper is to extend previous numerical work on rectangular cavities [22] to the case of a cylindrical 
one. 1 Previously [22] (for related work also see [23]), we considered the fundamental TM (TMm) and second 
fundamental TE (TE m ) mode for a semiconductor diaphragm (SCD) irradiated by a pulsed laser in a rectangular 
cavity. However, there are various reasons for considering a cylindrical cavity: the tuning of the standing wave 
frequencies depends only on the radius, R and length, L z (rather than L x , L y , L z ); it is easier to irradiate the SCD 
uniformly; and easier to construct a higher finesse cavity. 

Furthermore, in some proposed experiments, a Rydberg atom beam (which can be used to detect individual photons 
[24-26]) may lead more favorably to considering TM modes. 2 Thus, in this work we focus on the lowest excited (second 
fundamental) TMon mode (for a cylinder of length L z = 100 mm and radius R = 25 mm) that has a resonant pulse 
duration of T « 103 ps for the TM n mode with frequency /on = 4.83 GHz (cf. the TE m mode with T w 131 ps 
and /in = 3.83 GHz). 3 However, as we discuss in Section IV, TM modes are susceptible to polarization losses in the 
SCD coming from the backreaction of DCE created photons with the plasma sheet (this stems from the fact that the 
conductivity is related to the imaginary part of the dielectric function). Even so, the results presented here give upper 
bounds on the number of photons created, particularly given the fact that (see Section IV) decreased laser powers 
(implying less dissipation) still lead to significant particle creation. 

The outline of this paper is as follows. In the next section (Section II) we explain the plasma sheet model (Section 

II A), the necessary boundary conditions (Section II B) and the form of the eigenfunctions (Section II C). In Section 

III we discuss how to find the number of particles created using the Bogolyubov method (Section III A) and then 
show how the TMoio fundamental mode has no contribution to the photon creation rate (Section IIIB). In Section IV 
we argue on the dissipative nature of TM modes, based on polarization effects. Finally, analysis & discussion is given 
in Section V. In Appendix A the Hertz potentials approach to separate Maxwell's equations, including polarization, 
is discussed. 



1 An experiment has already been built at Padau University [15-17] using an SCD fixed to the wall of a rectangular cavity. 

2 Private communication with Prof. S. Matsuki. 

3 The TMqio case does not contribute to photon creation, see Section IIIB. 
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II. THE MODEL 

A. Plasma sheet model 

A microscopic model would be a more realistic way to discuss the coupling of DCE photons to the background SCD 
in a consistent way particularly to include dissipative effects (for particle creation in a crystal, based on a microscopic 
model without dissipation, see [6]). It is also interesting to note that working with a dielectric of finite thickness 
(resulting in more complicated Bessel functions, e.g., see [27, 28]), does not result in the same wave equation/junction 
conditions in the limit of an infinitely thin dielectric for TM modes [29]. 4 However, these issues are somewhat out 
of the scope of the current work and, for simplicity, we discuss the interaction of a plasma sheet irradiated by a 
laser in the hydrodynamic approximation, because it leads to semi-analytic expressions for the mode functions and 
eigenvalues. 

The Hamiltonian for a surface plasma of electrons of charge e and "effective mass" m* on a background electro- 
magnetic field is (using the minimal substitution) 

H = ^^[E-D + B-H] + J# x {-1—. (p 5 -en s A||) 2 + en SJ 4 ^(x-x s ) (1) 

where the canonical momentum is 

£ = (ps-en s (t)A ll )/m*n s (t) , (2) 

n s (t) is the "time dependent" surface charge density and Ay is the vector potential. The Hamiltonian, which is subject 
to constraints, implies that = 0, see Ref. [30], and thus, the electron momentum is related to the tangential vector 
potential by 

i = -eA ll /m% (3) 

which implies that the surface current density is 

K = en s (i)4 = -^#A|| . (4) 

Using surface continuity [31]: 

<7 + n-[VxnxK]=0 (5) 

where n is the unit normal (this is actually the charge conservation law: a + V • K = 0) along with the Lorenz gauge 
condition: 

d t A + V • A = (6) 

we find 

. = _e^ v ^ ^ a= ^M Ao , (7) 

m* m* m* 

where Aq is the scalar potential. Thus, in the plasma model we see that the surface charge density depends on the 
number of charge carries n s which can be made to vary in time by using a pulsed laser, with time profile n s (t), 
see Figure 1 (inset). To mimic current [15] and proposed experiments as closely as possible we model n s {t) by two 
Gaussian profiles t e and r joined smoothly to a plateau of length t (all of picosecond order), see Figure 1. 

B. Boundary conditions 

The boundary conditions for a charged plasma interface were derived very concisely in the work of Namias [31] 
(actually for charged moving interfaces), where for completeness we include the case where the interface is moving 
(v^O): 

(D 2 -Di)-n = c7 (B a -Bi)-n = 

n X (H 2 - Hi) - v • n(D 2 - D x ) = K n x (E 2 - E x ) - v • n(B 2 - B x ) = . 

(8) 



4 We have verified that TE modes do agree with the 'plasma sheet' model in the case of an infinitely thin dielectric [29]. 
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Here n is the unit normal pointing from a given region / into another region II, and t is any unit vector tangental to 
the surface. Although it might be interesting to consider how the mechanical vibrations of a two-dimensional electron 
layer affects photon creation, we set v = in the following (however see the discussion in Section IIIB). 

Substituting the relations for E, B, cf. equation (A6), into the above boundary conditions, with v = 0, we find the 
following continuity and jump conditions, e.g., see [22, 32]: 

disc y\ z=d = , disc d z $\ z=d = (9) 

disc d z ^\ z=d = ti^4^V(d) , disc $(d) = -AtS^0,3|,=d ■ (10) 

These equations can now be used to solve for the eigenvalues (see Section II C, below). This work focuses on a 
cylindrical cavity of radius R = 25 mm and length L z — 100 mm, see Figure 1. 



C. Wavefunctions and eigenfrequencies 

From the continuity and the jump conditions given above we have the following solutions for the wavefunction (for 
TE modes): 

,(TE) 



A m J \ sin (k mz z)v nm (x.x) 0<z<d 



where for a cylindrical section we have [33, 34]: 

tw(xx) = ~^=—, — A m g , (12) 

^ RJ n {y nm )^l-n 2 /yl m V RJ 

where y nm is the mth positive root of J' n {y) — 0, for TE modes. Then, due to symmetry the eigenvalue relation 
depends only on the z direction and reduces to the result we previously found [22]: 

sm{k mz L z ) _ e 2 n s (t) 

T 1,2 ™* V-L^J 



(fcmJ Tl sin(fc mz [L z - d]) sin(fc m ^d) k^ 



where the ± signs refer to TE and TM modes respectively (for TE modes drop the l/kj_ factor). 

In the above the wavefunction for TM modes is obtained by replacing sin — > cos and v nm (x±) — > r nm (xjj, where 

^nm(^l) ~~f^ 7 7 rd n |x nm — je ^ (14) 



/n RJ n+ i(x nm ) V R 

and a;„ m is the mth root of J n (x) = 0. 

For comparison with the rectangular case [22] we present some representative examples for the numerical solution 
of Equation (13) above, for the case of r) = 1/2 (the SCD placed at the midpoint) in Figure 2, for a driving period 
of 0(100) ps. The figure shows k n (t) for TE & TM modes, where as found in [22], TE modes up-shift in frequency, 
while TM modes down-shift. In contrast to the TE case, one important feature (also found for rectangular cavities) 
is that decreased laser powers still lead to large frequency shifts for TM modes (cf. the TE low power case in Figure 
2); this may be of relevance to dissipation, see Section IV. 

It may also be worth mentioning that TM modes depend explicitly on the transverse eigenfrequencies, kj_, (see 
Equation 13) and thus TM modes are influenced directly by the topology of the transverse section. 5 The physical 
reason for this is that the electric Hertz vector, II e , responsible for TM modes, leads to transverse magnetic Hm and 
perpendicular Ej^ fields, which induces electron motion parallel to the SCD and hence depends on the transverse 
dimension. TE modes on the other hand arise from the magnetic Hertz vector, II m where in this case the electric 
and magnetic fields are transverse and perpendicular, respectively: En, Hi. 



5 For actuation in time of a mirror or semiconductor in the longitudinal direction the cylindrical and rectangular mode functions k n (t) 
are identical for TE modes. However, the resonant frequencies and driving periods (as well as u) n (t): Equations (20, 27)) are different 
for a rectangular and cylindrical cavity. 
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FIG. 2. (Color online) Frequency variation for various k n for fundamental TEn„ modes and for second fundamental TMo\ n 
modes, with two different laser powers, in dimensionless units V mzx L z = 5000 (on left) and V mz ,^L z = 1 (on right). Numerically, 
for a given power (VmnxLz) and driving period T we carefully choose the profile n s (t) oc ti + to + i e , with Gaussians: ti = 
exp(— t 2 /2a'f ), i = e,r (denoting excitation and relaxation, respectively) by varying a e such that there is a smooth transition 
from relaxation to excitation when the pulse repeats. The profile is then substituted into Equation (13), which is solved for 
numerically, resulting in plots such as that above. 



III. PARTICLE CREATION WITHOUT LOSSES 



Before we begin discussing how to evaluate photon creation rate via the Bogolyubov method we would like to 
mention that most of the numerics in this work assumes a pulse train of order 10 pulses per train, which is ~ 1000 
ps for a driving period of order 100 ps. One possible benefit of shorter pulses per train would be that weaker laser 
powers are needed. However, recently, it was argued in [21] that our previous results for a rectangular cavity [22] 
require a larger number of pulses per train to obtain a photon creation rate of about 5/sec for TM modes, comparable 
to results for the TE mode with moving mirrors in a rectangular cavity [35]. Thus, in Figure 3 (see Section III A 
below) we also ran some simulations with ~ 30 pulses/train and found a photon creation rate of about 5/sec for the 
TMon mode. However, longer pulse trains may imply that dissipative effects start damping the photon creation rate 
(see discussion in Section IV) 6 . 



A phcnomcnological time dependent damping term can be included using the Heisenberg-Langevin approach [36], also see [28, 37]. 
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A. Instantaneous mode functions and the Bogolyubov method 

The quantum field operator expansion 

^( r ' *) = X] [ a m^™( r ' *) + aLV'mO: *)] (15) 

111 

of the Hertz scalars with instantaneous basis ansatz during irradiation is [4]: 

C ut (r, t)=J2 P^m(r, t), t > , (16) 

m 

where before irradiation t < (for TE modes) we have the standard stationary time dependence: 

^S(r,*) = -=J — sin — ^ «„ ro (x ± ) (17) 



(for TM modes replace sin — > cos, v„ m (xj_) — > r nm (xi) and (tp, VP) — > (0, $)). When Equation (16) is substituted 
into the equations of motion (on either side of, but not including the SCD located at z = d): 

vi# + ^*-a?tf = o , \/z^d, (is) 

(for TM replace ^ — > <I>) we find the following set of coupled second order differential equations:' 

oo oo 

P« + u? n (<)PW = -J2 [2M m „Pi s) + M mn P^ + M n iM me P^} (19) 

m I 

and for a cylindrical cavity (for TE modes) we have: 

jt)=c*[(^y+ki z (t)] (20) 



,2 



where y„ m is the n-th root of the Bessel equation J' m (x) = [33] (for TM modes see Section IIIB). Note the coupling 
matrix is defined by [4] 



Then given the scalar product: 



(0,v) = -i/ d A x(<t>r - <t>r) , (22) 

./ cavity 



the Bogolyubov coefficients are defined as: 



= (C\ C) , A™ = -(VC ut , t€T) (23) 



where in terms of the "instantaneous" mode functions it is possible to show that [8, 35] 



Urn „f„1 . / 1 



+ £ M, ro P, (n) ] , (24) 



where a mn is obtained by complex conjugation. The number of photons in a given mode (for an initial vacuum state) 
is then given by [38]: 

N m (t) = £ \(3 mn \ 2 (25) 



In the following we drop the x, y dependence of the mode functions and write the eigenfrequencies solely in terms of index m z = m. 
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Pulse duration [ps] Pulse duration [ps] 

FIG. 3. (Color online) Particle creation rate for 3000 ps (about 30 pulses) for the lowest TMoi n modes n = 1,2,3,4 with high 
& low laser powers, respectively: VL Z — 5000, 1. The number of mode couplings was cut off at ^ ma x = 71 and we found no 
discernible difference by increasing the number to £ ma x = 81. The driving period for both cases was set to T — 113 ps. Insets; 
Plot of unitarity bound, Equation (26), for the same parameters. 



where as a representative example, in Figure 3, we have plotted the lowest DCE created modes for the TM case with 
4 iax = 71 and for about 30 laser pulses (we checked convergence by going up to £ max = 81). The plot shows that in a 
given mode we can easily obtain a large number of photons ~ 5/sec and also shows that even modes make a greater 
mode contribution to particle creation. 

Interestingly, we see that on time scales of the order of 15 ~ 30 pulses the TM012 mode is more greatly enhanced 
for high laser powers. However, in this work we wish to focus on pulses trains of order 10 (up to about 1200 ps) and 
hence the TMrjn mode is typically more dominant (though not necessarily for decreased laser powers, see Figure 3 at 
right). 

In Figure 3 we also see that the total photon creation rate (sum of all mode contributions) increases by an order of 
magnitude. However, for times greater than ~ 1000 ps we find that the total number of modes needed to maintain 
convergence implies a much larger value of £ m ax ^ 0(100) and is out of the scope of the current work. 

In the following sections we shall solve for the amount of particle creation numerically with a cutoff at £ max ~ 50 
where we find the results do not change (for time durations up to about ~ 1200 ps). However, an independent check 
comes from the unitarity constraint [38]: 

e 

*-max 

^(|a m „| 2 - Wmn\ 2 ) = 1, (26) 

n 

which we have verified (e.g., see the insets in Figure 3 above, and also figures in References [22, 35, 39]). 8 



B. No particle creation for the fundamental TM mode 



We will now show that the fundamental TM io mode leads to no photon creation, because it is a zero mode and 
it also has a coupling matrix M mn = (this is converse to the case of moving mirrors [34, 40]. Thus, the second 
fundamental mode (TMon) is the lowest excited TM mode that produces DCE radiation in the plasma sheet model. 

To start consider the angular eigenfrequency for TM modes: 



».(*) 



R 



+ k 2 m (t) 



where x nm is the n-th root of the Bessel equation J rn (x) — [33]. The lowest eigenfrequency in the static case: 
with m,p = 0, 1, 2, . . . and n = 1, 2, 3, . . . becomes 

2 - 



' J oio 



(27) 



(28) 



8 More details of the numerical method can be found in [8, 22, 35, 39]. 
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and in the time dependent case instead of mode p 2 ir 2 /L 2 = 0, we have k 2 ^ (t) which also has a zero mode, k 2 ^ =0 (t) = 
0, which means that there is no parametric enhancement of the TMoio mode in Equation (19). Also note that by 
definition the coupling matrix M mn in Equation (21) is also zero and hence Equation (19) leads to no particle creation 
through multimodc coupling either. 

One might wonder, thus, how photon creation occurs at all for TMoio, even for the moving mirror case, given that 
we have such a zero mode? However, as discussed in [40] although the boundary conditions in Equation (8) (for p = 
and K = 0) lead to the standard Dirichlet condition for TE modes, a generalized Neumann boundary condition arises 
for TM modes, where for our Hertz potentials [34] we have: 

#0 = 0,Z)=0 (d + -dt)<5>(z = 0,L) = Q (29) 

c 

and we have assumed a perfect conductor with vanishing field in region II (the region external to the cavity for the 
moving wall case). 

This non-standard generalized Neumann boundary condition can lead to subtleties with quantization; however, by 
making a coordinate transformation we can work in a frame of reference where the time derivative vanishes. This 
leads to extra terms that appear in the coupled differential equation, Equation (19), e.g., see [40] and [34]. These are 
the terms responsible for DCE photon creation for TMoio modes for a moving boundary, but they are not present for 
the plasma sheet model (with v = 0), as we have discussed. 



IV. POLARIZATION LOSSES IN THE PLASMA SHEET MODEL 



It appears difficult to discuss dissipation in the plasma sheet model, because a (5-function profile has no intrinsic 
length scale in the longitudinal z-direction, where for example we expect temperature rises in the SCD to lead to 
dissipative effects. However, we can still make some general statements about dissipation based solely on Maxwell's 
equations: To consider dissipation we can either use a longitudinal z-dependent imaginary part in the dielectric 
function, £2(i), and set the electric polarization vector, P = 0; or conversely (see Appendix A) we can set £2 = and 
include a time dependent polarization P. 9 

In this work we will consider dissipation based on the assumption that losses arise due to electric polarization, 
P(t), 10 of pair created DCE photons, with back-reacted E(t). These photons are of GHz frequency (in the microwave 
regime) that implies that lot <C 1 where r is the typical relaxation time of a conductor. The laser field itself (which 
generates n s (t)) also leads to losses [21], but we will assume that most of the energy absorbed there is used to create 
the plasma window itself, i.e., to move the valence electrons into the conduction band. 

Based on the above assumptions we will now show that polarization losses from DCE photons only affect TM 
modes as follows: In Appendix A below Equations (A8) we are free to choose a gauge where all stream potentials 
are zero (except Q e ), see e.g., [42]. Thus, Maxwell's equations in Hertz form (Equations A8) imply that losses 
due to polarization only affect II e , namely TM modes 11 (assuming that (i(x,t) = po is a constant i.e., the induced 
magnetization is M = 0). 

In addition, for a cylinder the lowest frequency modes: (0,1,1) for TM and (1,1,1) for TE, give the dominant 
contribution to polarization, 12 because, quite generally, these are the modes with the greatest parametric enhancement 
due to resonance (& multimodc enhancement for TM modes, see Figure 5). Thus, in what follows we develop a 
simplified model of dissipation (through polarization effects) based on the Drude model. 

The electric polarization can be defined (at the linear level) as: 

P(t) = x(*)E(*) , (30) 
where in the Drude model the susceptibility in momentum space is, e.g., see [41]: 

xM = -— ,, , ■ (31) 

It is then possible to show that the real-time susceptibility (via an inverse Fourier transform) gives: 

*(*-*') = — e- (t - t,)/r , X(t-t')=0 for ?>t, (32) 



9 The two approaches here are equivalent, because the definition of the electric displacement is D = eE + P. 

10 Choosing P is equivalently to choosing — iojj 7^ 0, e.g., see [41]. 

11 Another way to understand this fact is that for transverse TE waves, Em and H^, have a lower order contribution to damping, e.g., see 
Chapter 8 of [33]. 

12 The TMqio mode is a zero mode and therefore not excited by DCE radiation, see Section IIIB. 
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where n v (t) is the volume charge density (related to the areal density n s {t) via the penetration depth, 8d- n v ~ 8dn s ), 
t is the relaxation (recombination) time and m* is the effective mass of the conduction electrons in the SCD. It 
may be worth mentioning that, because the polarization depends on the field strength we are left with a set of 
integro- differential equations in Equation (A8). 

To further simplify our discussion, we now assume that only the second fundamental TM mode (the lowest mode 
in our case) with won = 30.3 GHz gives the dominant contribution to polarization (denoted by ujq in what follows). 
A Fourier decomposition for a single mode implies (ignoring the fact that a bounded cavity would have sinusoidal 
modes) 

E(i) = E a e- lulot (33) 
and upon substituting this into the definition of causal polarization [41]: 



p oo 

/ dt" X {t")-E{t - 1") 
Jo 



P= / dt" X (t")~E(t-t") (34) 



along with Equation (32) we find 



P(t) = EoW2r / dt"e-*"/- e --o(t-t") ( 



E 5 d n s e 2 T 1 

7 — / V 1 

(I/t — Wq) 



(35) 



Hence, there are two ways to reduce losses due to polarization: One way is to decrease the laser power, because 
as discussed in [22], a laser power of 100 /iJ/pulse leads to a penetration depth, 84 ~ 50 /im; whereas for weaker 
laser powers, such as 0.01 ^iJ/pulse, we can reduce this depth by a factor of 100 ~ 1000. Interestingly, in the next 
section we will see that TM modes are excited for low laser powers (ignoring polarization), see Figs. 5 and 6. This 
compliments the fact that smaller values of 84 lead to less dissipation. 

The other way to reduce losses can be seen by taking the real and imaginary parts of 

(36) 



(1/t - iu ) 
in Equation (35). We see that the limit 

lo t < 1 (37) 

leads to Re[P] ^> Im[P]. Thus, for frequencies of Gigahertz order: loq ~ O(10) GHz, we can also reduce the amount 
of dissipation by reducing the recombination time in the SCD down to picosecond order: r ~ O(10) ps (this rather 
naive analysis leads us to conclude that relaxation times of nanosecond order are not sufficient, also see the discussion 
in [28]). Picosecond order relaxation times can; however, be achieved by an appropriate semiconductor doping, or by 
bombarding the SCD with Gold ions [43]. 

V. ANALYSIS & DISCUSSION 

Based on the assumptions just made above we can assume that under certain conditions, the Bogolyubov method 
(which does not include losses) will leads to results that give an upper bound on the particle creation rate. Thus, 
in Figures 4 & 5 we present results for the tuning dependence of the photon creation rate in a cylindrical cavity for 
TEni and TM n modes, respectively. The results are presented for the case where the SCD is placed at the midpoint 
for the lowest TE and TM modes in Figures 4 & 5, and as we found for a rectangular cavity, we find that multimode 
coupling enhances the TM contribution, while for the TE case, it diminishes (this is related to the different behavior 
of k n {t) for TE and TM modes, see Equation (13) and Figure 2). 

One of the most numerically intensive parts of this work is in the calculation of the position dependence, r\ = d/L z , 
of the photon creation rate. Except for the midpoint (rj — 1/2) the coupling matrix M mn (t) evaluates to thousands of 
lines of FORTRAN code and significantly slows down the numerics. 13 Furthermore, the procedure is hindered by the 



The current form of the equations are not parallelizable, and we solve for M mn (i) exactly for a given rj = d/L z 
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FIG. 4. (Color online) Tuning dependence for TEm modes, where the period of the pulse train varies from T = 119 ps up to 
T = 133 ps for a high laser power: VL Z = 5000 (resonant period is 131 ps). The particle creation rate is that at t ~ 1200 ps. 



FIG. 5. (Color online) Tuning dependence for TMoii modes; red circles and black sguares are for high powers with single 
and multimode coupling, respectively. Blue diamonds and green triangles are for low laser powers, with single and multimode 
couplings, respectively. The period of the pulse train varies from T = 102 ps up to T — 115 ps (resonant period is 103 ps). The 
particle creation rate is that at t ~ 1200 ps. 

fact that the coupling matrix detunes the resonant driving period of the laser pulse train and we have to make many 
runs at different driving periods (this is also an effect we expect from dissipation as well [36] ) . Hence in this work we 
focused on a given probe driving frequency of T = 113 ps, except for cases where we found no real enhancement. 1 
In Figure 6 we have plotted the position dependence for TMon modes, and like for a rectangular cavity we find that 
TM modes are enhanced even for decreased laser powers. They are also quite generally unaffected by the location of 
the SCD, as compared to TE modes. 

In summary, we have evaluated particle creation rate in a microwave cylindrical cavity, where a laser periodically 
irradiates a semiconductor diaphragm (SCD). We focused on the second fundamental TMon mode (where the time 
dependence of the SCD was modeled using the 'plasma sheet' model, ignoring dissipation). Importantly we showed 
that TMon modes are fairly insensitive to the position, 77, and also produce significant numbers of photons for 
decreased laser powers, as opposed to TE modes. Thus, because we have so far ignored dissipation, our numerical 
results should give upper bounds on the amount of photons created for TMqii modes. 

We also explained in Section III B, why the TM010 is not excited for a cylindrical cavity in the plasma sheet model. 
In the case of moving mirrors and, for example, working with an instantaneous basis, see Equation (16), we obtain 



For rj = 0.4 we found that T = 107 ps was a better driving period in Figure 6. 




103 105 107 109 Ii4 1% 
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Position dependence: n=d/L z 

FIG. 6. (Color online) Position dependence for TMqh modes (with multimode coupling); red solid and blue dashed lines line 
are for high (VL Z = 5000 ) and low (VL Z — 1) laser powers, respectively. The frequency of the pulse train was chosen to be 
T = 113 ps for all cases except for high power, n = 0.4 case where T — 107 ps (resonant period is 103 ps). The particle creation 
rate is that at t ~ 1200 ps. 



extra terms, besides the coupling matrix, M mn , and zero mode, cl> io: which are not zero for TM io- This is what 
leads to particle creation in this mode, e.g., see Equation (37) in Ref. [34]. However, in the case of an irradiated 
plasma sheet the boundary conditions are different and no extra coupling terms arise: the TMoio is not excited. 

In Section IV we discussed dissipation from the electric polarization of DCE pair created photons. An important 
point from this analysis shows (at least within the plasma sheet model) that only TM modes are affected, while TE 
modes are not. However, we also argued that decreased laser powers reduce polarization, which is encouraging given 
that decreased powers still lead to significant photon creation (for TM modes), see Figure 5 and 6. 

Of course there are limitations and we are currently developing a numerical method to include the effects of 
dissipation for TE and TM modes in a rectangular/cylindrical cavity [36]. This approach leads to a parametric 
equation with detuning very similar to that found via the Bogolyubov method with multimode coupling, see Equation 
(19), and may serendipitously result in TM modes enhancing in a way that leads to an asymptotic saturation of 
dissipation (see [44] for a discussion of this point for single mode coupling). Furthermore, as we mentioned earlier 
there are limitations to the validity of the 'plasma sheet' model and indeed it remains to be seen if the behavior of 
k n (t) for decreased laser powers continues in some kind of microscopic model of a plasma sheet (or by modeling the 
SCD as a thick dielectric [27, 28]). These issues as well as higher mode effects are left for future work. 

Hopefully, the work presented here leads support to the idea that TE and, in particular, TM modes are well suited 
to the detection of DCE radiation in a cylindrical cavity (even with possible losses) . We hope these results may be of 
use for current and proposed experiments to detect DCE radiation in microwave (centimeter-sized) cavities. 
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Appendix A: Hertz potentials approach 

For completeness we describe, following Refs. [22, 34], how to use Hertz vectors to define a set of potentials, e.g. see 
[33, 42], which conveniently separate Maxwell's equations into TE and TM equations of motion. This allows one to 
essentially work with two scalar field equations (with different boundary conditions). In the following we shall review 
the discussion given for example in [33, 34, 42]. 
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In the Lorenz gauge: 

[ieA + V • A = , (Al) 
where A is the scalar potential Maxwell's equations become: 

^5#- V ^o = V-V-Po (A2) 
ot e e 

^-V 2 A=-J + ^+VxM„ (A3) 
at (i ot 

where the permanent polarization and magnetization (Po, Mo) are introduced to motivate the form of the potential. 
By defining two Hertz vector potentials II e and U m as (with fiQ = 1) 

A = --V-U e A = /i^ + Vxn m (A4) 

then is it straightforward to show that Equations (A2, A3), above, automatically satisfy the Lorenz gauge condition, 
Equation (Al), as can be verified. Then from the definition of the electromagnetic field in terms of 

B = V x A , E = -d t A - VA Q (A5) 

the electric field and magnetic displacement can be written in terms of Hertz vectors as 

E = -v ( v.n t) -^-vx^, 
on 

B = /iVx + V x (V x n m ) . (A6) 
The separation is effected by introducing the following stream potentials [42] 

p=-V-Q e , J = Q e + -VxQ m (A7) 

M 

with a similar result for the magnetic stream potentials R e and R m , which have zero magnetic charge and current. 
Using these definitions, Maxwell's equations separate into [42]: 

Aie ri e - v 2 n e = (p + Q e + R e ) 

HE ri m - V 2 n m = (M + Q m + R m ) . (A8) 

From the theory of gauge transformations of the third kind [42] 15 it is always possible to choose a gauge where the 
stream potentials are Q m = R e = R m = and Q e = J dtJ and in the 'plasma sheet' model J = (no bulk charges, 
only surface charges, K). 

From the symmetry involved it is convenient to define the following Hertz potentials: 

n e = $ z n m = * z , (A9) 

where z is a unit vector in the longitudinal direction z and $ and ^ correspond to TM and TE modes respectively. 
It is then easy to show that 

A Q = -^d z <Z>, A = y *x-0 x *y + M t $z . (A10) 

The potentials ^ and $ represent TE and TM modes, respectively and similarly the E and B fields become (using 
equation (A6)): 

E = (^0,0,* - d y d t ^)± + {- e d y d z <S> + d x d t V)y - -V\cj> z 

B = {»d y d t $ + a^tf)x+ (-iid x dt$ + d y d z ^)y - Vi^z (All) 



Nisbet in [42] defines gauge transformations on 2-form fields as transformations of the third kind. 
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where for cylindrical coordinates (p, 9, z) we have the following transverse Laplacian 




) 



+ 



(A12) 



The great utility of the Hertz potentials approach is that the separation leads to a set of two scalar field equations 
(one for TE and the other for TM) where the longitudinal symmetry (the axis of a cylinder or rectangle) is decoupled. 
Thus, considering 1 + 1 dimensions or 3 + 1 dimensions does not really complicate the problem (see [35] for TE modes 
in 3 + 1 dimensions) . 
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